Strong mass effect on ion beam mixing in metal bilayers 
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Molecular dynamics simulations have been used to study the mechanism of ion beam mixing in 
metal bilayers. We are able to explain the ion induced low-temperature phase stability and melting 
behavior of bilayers using only a simple ballistic picture up to 10 keV ion energies. The atomic mass 
ratio of the overlayer and the substrate constituents seems to be a key quantity in understanding 
atomic mixing. The critical bilayer mass ratio of S < 0.33 is required for the occurrence of a thermal 
spike (local melting) with a lifetime of r > 0.3 ps at low-energy ion irradiation (1 keV) due to a 
ballistic mechanism. The existing experimental data follow the same trend as the simulated values. 

PACS numbers: 79.20.Rf, 61.80.Az 61.80.Jh 61.82.Bg 
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The mechanism of atomic transport in solids is a fun- 
damental subject of material science [1]. The atomic re- 
locations, such as atomic mixing at interfaces under the 
effect of ion irradiation has also been an important topic 
in the last decades [2,3]. In these phenomena the atomic 
systems are driven far from equilibrium and the roughen- 
ing of the interface acts in competition with the restoring 
forces of the relaxation process. In certain interfacial sys- 
tems ion beam mixing (IM) or the elevated temperature 
seems to randomize the initially sharp interface leading 
to amorphization [4,5], to crystal-glass [3] or liquid to 
glass transition [6] . In many bilayers no considerable IM 
occurs [7,8]. As an explanation for the ion induced phase 
stability (instability) of interfaces, the thermal spike (TS) 
model is widely accepted in the last two decades [8,9] 
which predicts the dependence of the mixing efficiency 
on thermodynamic quantities such as heat of mixing. Re- 
cently though we observed the presence of TS but failed 
to find any effect of heat of mixing on IM [10,11]. Thus 
accepting the importance of the TS we propose to under- 
stand its affect on IM as a purely ballistic phenomenon 
[11]. 

In the present Letter we would like to show that 
the atomic mass ratio has a dramatic effect on IM in 
various metal bilayers. At 1 keV Ar + ion bombard- 
ment no real TS occurs above the atomic mass ratio of 
S — m over i/mbuik > 0.33 while considerable IM occurs if 
S < 0.33, where m over i and rribuik are the atomic masses 
in the overlayer and in the substrate (bulk). We will show 
that the measured mixing efficiencies follow the same 
trend as the simulated IM. We propose to understand 
IM as a ballistic process and we find that the chemical 
interdiffusion model [12] might be inappropriate. 

Classical constant volume molecular dynamics simula- 
tions were used to simulate the ion-solid interaction using 
the PARC AS code [13]. The computer animations can be 
seen in our web page [14]. Here we only shortly summa- 
rize the most important aspects. Its detailed description 
is given in [13] and details specific to the current systems 
in recent communications [10,15]. We irradiate a series of 
bilayers (Al/Pt, Ti/Pt, Al/Ag, Cu/Pt, Cu/Au, Al/Cu, 



Ni/Ag, Ag/Au, Ti/Co, Ag/Ni, Cu/Ni) with 1 keV Ar+ 
ions. In certain bilayers (Al/Pt, Ti/Pt, Al/Cu, Ni/Ag, 
Ti/Co) we increase the ion energy up to 10 keV in order 
to show that the mass effect persists up to higher ion 
energies. The initial velocity direction of the impacting 
atom was 7 degrees with respect to the surface of the 
crystal (grazing angle of incidence). To obtain a repre- 
sentative statistics, the impact position of the incoming 
ion is varied randomly within a 5 x 5 A 2 area. 

We used a tight binding many body potential given 
by Cleri and Rosato (CR) [16], to describe interatomic 
interactions. We have chosen those bilayers for which 
atomic potentials are available. This type of a potential 
gives a very good description of lattice vacancies, includ- 
ing migration properties and a reasonable description of 
solid surfaces and melting [16]. The bilayer AB cross po- 
tential is composed of the CR potentials of the A and 
B elements. Taking the simple average of the elemental 
potentials the heat of mixing is set to AH m w 0. This 
can be done since we have shown recently, that we found 
no considerable dependence of IM on AH m in Ti/Pt [10] 
up to 7 keV ion energy [11]. The choice of AH m « 
is also useful from that point of view that we would like 
to show that even in the case of AH m rs (the lack 
of AB chemical bonding) considerable IM occurs in cer- 
tain bilayers. In other bilayers the effect of AH m on IM 
is also very weak if any [11]. The reliability of the AB 
crosspotential as an average of the elemental potentials 
is tested in the case of CuAu, where an optimized CR 
potential is available for the alloy phase [17]. We find no 
considerable changes in the physical properties or in the 
equilibrium structure. Although this does not prove the 
reliability of other crosspotentials, it suggests that the 
average potentials are very close to the optimized one. 

The construction of the interface systems is given else- 
where [11]. We only shortly summarize that the inter- 
faces have (111) orientation and the close packed direc- 
tions are parallel. The thickness of the upper layer is 4-8 
monolayer (ML), while the bulk is constructed from 36 
MLs in those samples which subjected to 1 keV irradia- 
tion. These samples includes roughly 45000 atoms. At 
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FIG. 1. The number of intermixed atoms (N m i x , interfacial 
mixing) at various atomic mass ratio (<5) in various metal bi- 
layers obtained at 1 keV Ar + irradiation. Note the threshold 
8 « 0.33 value below which the bilayers exhibit an ion-beam 
mixing "catastrophe" and above which only cascade mixing 
occurs. The data points correspond from left to right to the 
bilayers Al/Pt (0.14), Al/Ag (0.24), Ti/Pt ( 0.25), Cu/Au 
(0.32), Cu/Pt (0.33), Al/Cu (0.42), Ni/Ag (0.54), Au/Ag 
(1.83), Ti/Co (0.81), Cu/Ni (1.08) and Ag/Ni (1.82), respec- 
tively (the atomic mass ratios are given in the paranthesis). In 
those case where 5 > 1 we use 1/5, because we find that N m i x 
in AB and BA systems are nearly equal. The error bars de- 
note standard deviations. Inset: The number of mixed atoms 
is also shown for higher energies for Al/Pt (6 keV), Ti/Pt (8 
keV), Al/Cu (10 keV), Ni/Ag (9 keV) and for Ti/Co (10 keV) 
as a function of the mass ratio. 



higher irradiation energy we put a thicker overlayer with 
8-16 MLs and a substrate with 90 MLs (550000 atoms). 
The interfacial systems are created as follows: the over- 
layer is put by hand on the (111) substrate (bulk) and 
various structures are probed and are put together ran- 
domly. Finally that one is selected which has the smallest 
misfit strain prior to the relaxation run. The remaining 
misfit is properly minimized during the relaxation pro- 
cess so that the overlayer and the substrate layers keep 
their original crystal structure and we get an atomically 
sharp interface. 

We carry out liquid and hot atom analysis in order to 
study the role of highly mobile particles in IM [11]. To 
do so we need to calculate the time dependent temper- 
ature of the ith high energy particle {Ti t i oca i{t)) which 
can be given as follows: ^rriivf(t) — ^kTij ocal (t) where 
rrii and Vi are the atomic mass and the velocity of the 
ith high energy particle and k is the Boltzmann con- 
stant. We allow some temperature fluctuation for the 
liquid atoms above T m (the melting point), hence for a 
hot atom and for recoils we use the arbitrary definition 
T m + 1000 < Tiorai [11]. Hot atoms are those energetic 
particles which have a medium long mean path (~ 10 A) 
with an average kinetic energy of 0.1-1 eV. Liquid atoms 
have shorter, recoils have much longer mean free path. A 
liquid atom can be considered as a nearly standing atom 
during the lifetime of the TS and only vibrationally ex- 
cited around its equilibrium position while a hot atom 
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FIG. 2. The low-temperature experimental mixing effi- 
ciency (A 5 /eV) as a function of the atomic mass ratio (8). 
The values are taken from refs. [7,8]. In those cases where 
8 > 1 we use instead 1/8 because we find that there is no seri- 
ous difference in mixing between AB and BA bilayers. Inset: 
The calculated (simulated) ^im (A 5 /eV) vs. 8. 



jumps to another position with a velocity of Vi w 1 A/ps 
[11]. The hot atoms are present till the end of the TS 
hence cascade and TS mixing coexists in certain interfa- 
cial systems. The more precise definition of a hot atom 
is given elsewhere [11]. 

We calculate the number of mixed atoms (N mix ) and 

the simulated mixing efficiency = g^gj; , where 

(R 2 ), no and Eu N , are the calculated mean square atomic 
displacement through the interface per atom, the atomic 
density in the upper layer and the deposited nuclear en- 
ergy. We exclude from (R 2 ) atomic displacements which 
do not lead to broadening at the interface (self-atomic 
mixing) because the experimental £jm is calculated from 
broadening at the interface [8,9]. Further calculational 
details are given in ref. [18]. 

We investigate the influence of the atomic mass ratio 
on ion beam mixing. We ion bombarded various bilayer 
systems with different atomic mass ratios and find that 
below a threshold ratio (5 < 0.33) the magnitude of in- 
termixing (N mix ) is enhanced abruptly. The results are 
summarized in FIGs (l)-(2). In FIG (1) we plot the sim- 
ulated N m i x vs. 5 and the experimental £/m vs. <5 in 
FIG (2) collected from refs. [7,8]. In the inset FIG (2) we 
also give the simulated £jm vs. <5 for bilayers for which 
high energy (up to 10 keV) simulations are available and 
also for the other bilayers for which is calculated at 
1 keV ion energy. In FIG (2) wee see that the increase 
in £im occurs between <5 m 1/3 and 1/2 which is rather 
similar to the inset FIG (1). In FIG (2) the experi- 
mental mixing efficiency values are obtained as follows 
[7]: (£/m ~ ir^)) where Dt is the diffusion length, <f> 
is the fluence and Fjj is the deposited energy at the in- 
terface). Although there is some scatter in the data at 
around S « 1/3, £jm increases heavily in accordance with 
the N mix values shown in FIG (1). Although the values 
in FIG (2) are obtained at high energies, however, they 
can be compared with our £|™ values since IM occurs 
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FIG. 3. The xz atomic positions of the recoils and the hot 
atoms at 1 keV ion energy collected up to 4 ps in Ti/Pt. The 
z-coordinate is the depth position. The dashed line denotes 
the interface. 



FIG. 4. The xz positions of the recoils and the hot atoms 
at 1 keV ion energy collected up to 0.3 ps in Ti/Co. The 
z-coordinate is the depth position. The dashed line denotes 
the interface. 



primarily in the subcascade region, where the energy is 
in the range we applied [18]. 

In particular in Al/Pt we find a strong amorphization 
and broadening at the interface which is in accordance 
with measurements [4,5,11,19]. In Cu/Pt and in Cu/Au 
the phase stability of the interface is also very weak in ac- 
cordance with ion irradiation experiments [19]. In Cu/Au 
1 MeV ion bombardment results in strong broadening 
[8,19]. In bilayers Al/Ag and in Au/Ag we find a rel- 
atively weak interfacial mixing (N m i X < 30), however, 
is large (£f£? w 300 ±120 in Al/Ag, £f£? « 115 ±45 
in Au/Ag). A large measured value of approx265 
is found in Au/Ag [7] and this data point is plotted at 
1/5 = 0.55 in FIG (2). Interestingly the atomic mobility 
is relatively large in these systems while the gross num- 
ber of mixed atoms (N m i x ) is small. We attribute this 
" anomalous" behavior of Au/Ag to the tendency of crater 
formation in Au [9]. It has been shown, recently, that 
crater formation enhances mass transport hence atomic 
mobility between the interface and the free surface [15]. 

In FIG (2) the experimental mixing efficiency values 
[7,8] (£/m ~ ^j)> where Dt is the diffusion length, $ is 
the fluence and Fd is the deposited energy at the inter- 
face) are plotted against the atomic mass ratio. Although 
there is some scatter in the data at around S « 1/3, 
£tm increases heavily in accordance with the N mix val- 
ues shown in FIG (1). Although the values in FIG (2) 
are obtained at high energies, however, they can be com- 
pared with our ^f^p values since IM occurs primarily in 
the subcascade region, where the energy is in the range 
we applied [18]. 

For bilayers with 6 > 0.33 starting with Al/Cu (5 = 
0.42) we got a weak interfacial mixing in accordance with 
the MD [20,21] and experimental studies [7]. The strong 
mass effect can be understood as the ballistic mechanism 
plays an essential role in IM and may stem from an in- 
creased backscattering of light overlayer atoms from the 
heavy substrate [11]. Moreover the density of the colli- 
sional cascades depends on the bilayer mass ratio [11]. 



It should also be stressed that in pure elements or even 
in metal alloys we find a much weaker IM and shorter 
TS [14,15] due to the strong (re)ordering forces. In order 
to elucidate the relative roles of ballistics and TS in the 
interfacial mixing we examined the prototypical bilayer 
Ni/Ag (S « 0.54) and found a very weak interfacial mix- 
ing at 1 keV bombardment and no occurrence of a real 
TS period. The energetic atoms (recoils, hot atoms) dis- 
sappear at less then 0.3 ps, which is typically the end of 
the cascade period. This system is relatively well studied 
theoretically [20] at 10 keV ion energy and being a typical 
example of a segregating system which has a high posi- 
tive heat of mixing in the liquid [20] . At 9 keV we find a 
real thermal spike which persists up to 3 ps. The number 
of mixed atoms (N m i X ) is, however, much smaller then in 
Al/Pt or in Ti/Pt (inset FIG 1)). Therefore, although 
there is an increase in N mix , the trend remains the same: 
the effect of the mass ratio is also robust at higher ener- 
gies. We attribute, however, the weak IM in Ni/Ag not 
to the positive heat of mixing, but to the weak mass ef- 
fect in this system (0.33 < S w 0.54). In Ti/Co we find 
no TS up to 10 keV as well as IM is very weak in this sys- 
tem (5 w 0.82). It should be emphasized that the similar 
situation is true for all the bilayer samples which have 
S > 0.33 (FIG (1)). This observation clearly indicates a 
robust mass effect when 6 < 0.33, hence a kinematic and 
ballistic picture seems to be sufficient for describing IM 
in bilayer systems. In those systems, where TS does not 
occur, cascade mixing is the only IM effect. 

The cascade events during the ballistic period is shown 
in FIGs (3)- (4). In the case of Ti/Pt we get a dense col- 
lision cascade, the recoils (hot atoms) are concentrated 
within a smaller region due to couple of reflections while 
in Ti/Co the high energy particles are scattered in a 
larger volume hence the deposited energy spreads over 
a larger irradiated region. Therefore if the mass ratio 
S > 0.33, the cooling of the cascade is ultrafast due to 
the low concentration of the recoils (hot atoms). Indeed, 
we found that the average atomic concentration of the 
hot atoms in the irradiated zone (V w 1000 A ) is around 
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10 22 atom/cm 3 for S > 0.33 and 4 x 10 22 atom/cm 3 for 
S < 0.33. These values should also be compared with the 
average atomic concentration of 5 — 8 x 10 22 atom/cm 3 
in metals. If the mass ratio drops below 0.33, the hot 
particles are still present in the TS. One can see that 
for S < 0.33 the peak hot atom concentration is close to 
the atomic concentration. In these systems we no longer 
have a simple liquid ensemble, it is rather a superheated 
system [11]. The mass effect is even more pronounced 
at higher energies. The corresponding plots of the colli- 
sional cascades in Al/Pt at 6 keV ion energy can be seen 
in a web page [22] . We have also shown that the density 
of the cascade is strongly affected by <5 [11,14]. 

Qualitatively we explain the observed strong mass ef- 
fect as follows: In elastic collisions of the recoils (ener- 
getic light particles from the overlayer) with the heavier 
substrate atoms the kinetic energy of the moving atoms 
is partly transferred to the heavier atoms, which , how- 
ever, might not be kicked out of their positions because 
of the large mass difference. The colliding heavy part- 
ner of the recoil becomes vibrationally excited, which 
means that its rms amplitude of thermal vibrations be- 
comes equal to about 50 % of the interatomic distance. 
That is basically the Lindeman's criterion for lattice in- 
stability during melting: a crystal melts when the rms 
thermal displacement of atoms from their equilibrium 
positions become large enough to invade their nearest- 
neighbor spaces [3]. For such thermal displacements, the 
thermal expansion would far exceed the critical value for 
shear instability (the Born criterion, the lost of at least 
one of the shear moduli) leading to mechanical instability 
[3]. The neighborhood of this hot heavy atom is heated 
up and local melting occurs (TS). When the mass of a 
recoil and the colliding partner is comparable the target 
atom might be displaced from its original position leav- 
ing a vacant site. In this case the slowing down of the 
recoil in the bulk does not result in local melting because 
the kinetic energy of the recoil spreads over a too large 
volume. 

The important question remained to be answered what 
is the reason of the critical mass ratio of <5 ~ 0.33? We do 
believe that the threshold value is due to the emergence of 
a strong backscattering of the recoiling light atoms at the 
interface. Below this value hence the energy deposition 
becomes extremely effective at the interface through en- 
ergy transfer to the standing heavy atoms. The backscat- 
tering effect at the interface results in the confinement of 
the light recoils in the overlayer which leads to superheat- 
ing. This is the primary reason of the high concentration 
of hot atoms in these bilayers (see FIGs (3)). The interra- 
cial backscattering phenomenon can be attributed partly 
to the mass difference and also to other effects such as 
the difference in the cohesive energies in the substrate 
and in the overlayer [11]. 

In summary, we have shown that intermixing in metal 
bilayers strongly depends on the relative masses of the 
constituents under the effect of ion irradiation. There 
exists a threshold mass ratio value below which the in- 



terface system is unstable against ion bombardment. We 
propose to understand ion beam mixing as a ballistic 
process. The observed strong mass effect in heterophases 
might be an important topic in preparation of thin films 
and multilayers especially with great technological im- 
portance. 
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